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Abstract. We consider the Unipolar Inductor Model (Goldreich & Lynden-Bell 1969) applied to Double Degenerate 
Binaries (DDEs) with ultrashort periods (Wu et al. 2002). 

In this model a magnetized primary white dwarf has a slight asynchronism between its spin and orbital motion, 
so that the (non-magnetic) secondary experiences a motional electric field when moving through the primary field 
lines. This induces a current fiow between the two stars and provides an electric spin-orbit coupling mechanism 
for the primary. 

We study the combined effect of Gravitational Wave emission and electric spin-orbit coupling on the evolution 
of the primary degree of asynchronism and the associated rate of electric current dissipation in such systems, 
assuming that the primary's spin is not affected by other mechanisms such as tidal interactions with the companion. 
In particular, we show that in ultrashort period binaries the emission of GW pumps energy in the electric circuit 
as to keep it steadily active. This happens despite the fact that spin-orbit coupling can rapidly synchronize 
the primary, because GW represent a slow desynchronizing mechanism that steadily substracts orbital angular 
momentum to the system. A slightly asynchronous steady-state is thus achieved, determined by the balance 
between these two competing effects. This can be shown to correspond to a condition where the total available 
electric energy is conserved, because of GW emission, while dissipation, synchronization and orbital shrinking 
continue. 
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1. Introduction 

The Unipolar Inductor Model (UIM from here on) has 
been originally proposed (Goldreich & Lynden-Bell 1969) 
for the Jupiter-Io system to explain the origin of bursts 
of decametric radiation received from the planet, whose 
occurrence and intensity were known to be strongly influ- 
enced by the orbital location of lo. 

The model relies on the fact that Jupiter's spin differs from 
the system orbital period, lo's spin being locked to the or- 
bit around the planet. Given the '--^ 10 G surface dipole 
magnetic field of Jupiter and the good electrical conduc- 
tivity expected from the satellite, the system behaves as a 
remarkably simple electric circuit, that can be sketched as 
follows. A good conductor (the satellite) has a trasverse 
motion with respect to the field lines of the planet and 
this induces an e.m.f. across it^. The e.m.f. can acceler- 
ate free charges that are present in the ambient medium. 
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^ It would be more appropriate to describe this in terms of 
the Lorenz force acting on the charge carriers within the con- 
ductor, but introducing the e.m.f. emphasises the analogy with 
an electrical circuit 



giving rise to a flow of currents between the two objects. 
Currents are confined to a thin sheet along the sides of the 
flux tube connecting Jupiter and lo; hence, the cross sec- 
tion of the current layers at Jupiter's surface has the form 
of two arc-shaped strips. Currents flow along one side of 
the tube towards Jupiter and then vertically through its 
upper atmosphere and ionosphere. Jupiter's electrical con- 
ductivity (ct) in the ionosphere is expected to be almost 
isotropical (Goldreich & Lynden-Bell 1969 and references 
therein): hence, currents reaching it can propagate trans- 
verse to field lines, closing the circuit and returning back to 
lo along the opposite side of the fiux tube. Charges along 
the fiux tube are accelerated to mildly relativistic ener- 
gies and lead to coherent cyclotron emission over a range 
of wavelenghts: this is the framework in which Jupiter's 
decametric radiation and its strong modulation by lo's po- 
sition can be explained. 

Among the main expectations of the model it is of inter- 
est here mentioning the effect of resistive dissipation of 
currents in the planetary atmosphere: this causes a local 
heating and an associated localized enhancement of the 
thermal emission. 

Several confirmations to the model have been obtained 
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over the years, the most spectacular being provided by the 
HST UV observations of Clarke et al. (1996). In particular, 
these revealed the emission from the resistive dissipation 
region on Jupiter's surface - lo's footprint - and its coro- 
tation with the satellite, which means that the emission 
region drifts on the planet's surface exactly as expected in 
the UIM. 

The same basic picture has been proposed by Li et al. 
(1998) for planetary companions to white dwarfs: if the 
UIM applied to such systems it would offer an alternative 
way of searching for extrasolar planets through the elec- 
tromagnetic emission associated to the electrical circuit. 
Wu et al. (2002) proposed a similar scenario for the case of 
two white dwarfs forming a close binary system. They con- 
sider a moderately magnetized primary, a non-magnetic 
secondary and a primary spin not perfectly synchronous 
with the orbital motion, while the secondary's is efficiently 
kept synchronous by tidal forces. The atmospheric electri- 
cal conductivity of white dwarfs is expected to be highly 
anisotropic, because of the WD magnetic field, up to at 
least ^ 1 Km depth (Li. Forrario and Wickramashigc 
1998), so currents are forced to follow field lines in this 
situation: they can only cross them and close the circuit 
when (T becomes isotropical and cross-field diffusion can 
proceed efficiently. 

This model was proposed to account for the observed 
properties of the two candidate ultrashort period binaries 
RX J0806-I-15 and RX J1914-h24. These have significant 
soft X-ray emission (~ 10^^ and 10^^ erg respectiv- 
ley) pulsed at periods ^321 and 569 s: these periodicities 
are interpreted as likely due to orbital modulations, al- 
though the subject is still debated (see Cropper et al. 2003 
for a review on this and several other aspects). Further, 
both sources have a measured orbital spin-up (Strohmayer 
2004, Hakala et al. 2004, Israel et al. 2004), as if driven by 
GW with no matter transfer: this would thus require an al- 
ternative source for the X-ray emission, UIM representing 
an interesting possibility. In fact, given the very short or- 
bital periods, possibly larger primary magnetic fields and 
the compactness of both components with respect to pre- 
vious versions of the UIM, such systems may dissipate 
energy at a significant rate even for slight degrees of asyn- 
chronism, as shown by Wu et al. (2002). These authors 
showed that currents are resistively dissipated essentially 
in the primary atmosphere, the associated heating causing 
the observed soft X-ray emission: its source is ultimately 
represented by the energy of the relative motion between 
the primary spin and the orbit. Further, the process re- 
distributes angular momentum - via the Lorentz torque 
on cross-field currents - between the primary spin and the 
orbit and thus acts to synchronize them. 
In this application of the UIM, synchronization timescales 
~ few 10^ yrs were obtained for both RX J1914-I-24 and 
RX J0806+15, very short compared to the orbital evolu- 
tionary timescales, respectively ~ 8 x 10^ yrs and 3x10^ 
yrs. Accordingly, a very low probability (significantly < 
1%) of detecting these systems during the asynchronous 
phase would be expected. This would in turn require 



a very large population of such systems in the Galaxy, 
much larger than predicted by population-synthesis mod- 
els, since two of them have been detected in the short- 
lived, active phase. 

Concerning this point, we focus in this work on a key 
aspect of the problem that has been overlooked in pre- 
vious works, namely that in the framework of the UIM 
applied to DDBs perfect synchronization between the pri- 
mary spin and the orbit is never reached: therefore, the 
current flow should not stop at any time. This happens 
even if the synchronization timescale (r^) is much shorter 
than the orbit evolutionary timescale (tq — lOo/i^o), be- 
cause of the continuous loss of orbital angular momentum 
caused by GW. As tidal synchronization is expected to be 
efficient only for the lighter star, GW drive the primary 
out of synchronism on the longer timescale Tq, thus con- 
tinuosly feeding energy to the electric circuit. 
Consider a perfectly synchronous system: the electric cir- 
cuit would be turned off while GW would still cause orbital 
spin-up and, thus, a desynchronization of the primary's 
spin, which would in turn switch the circuit on: the sys- 
tem should thus evolve, over the timescale r^, towards a 
slightly asynchronous steady-state that is determined by 
the balance between the fast, synchronizing mechanism 
(UI) and the slow, desynchronizing one (GW). 
In the following this issue will be addressed quantitatively: 
in particular we show that, in a binary system with UIM at 
work, the orbital period decreases because of GW emission 
and the primary spin is forced by the coupling to approach 
the orbital period, but perfect equality is not reached be- 
cause of the energy fed by GWs to the circuit. Further, 
we derive a definition of Tq, as a function of system pa- 
rameters, find conditions under which Tq <C Tq is verified 
and discuss the salient evolutionary features of the UIM 
implied by this condition. 

In a companion paper we apply these general considera- 
tions - and the related formulas that we derive here - to 
RX J0806-I-15 and RX J1914-I-24 and obtain interesting 
constraints on system parameters for the UIM to be ap- 
plicable to these sources. 

Detailed and systematic calculations and evolutionary im- 
plications are beyond the scope of these works and will be 
addressed elsewhere. 

2. Asynchronous Evolution in the Unipolar 
Inductor Model 

Following Wu et al. (2002), the primary asynchronism pa- 
rameter is defined as a = coi/uio, where wi is the primary 
spin and Uo the orbital motion. 

Given an asynchronous system with orbital separation a, 
the secondary star will be moving across the primary mag- 
netic field lines with the relative velocity v = a{ijJo — coi) = 
[GMi{l + g)]i/3 col^^{l - a), where G is the gravita- 
tional constant. Mi the primary mass, q = M2/M1 the 
system mass ratio. The electric field induced through 
the secondary is E = ^ ^ ^ , with an associated e.m.f. 
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$ = 2R2E, R2 being the secondary's radius. Because of 
the induced e.m.f. a current flows between the two compo- 
nent stars, whose resistive dissipation - which takes place 
essentially in the primary atmosphere - has two effects: 
first, it causes significant heating of the dissipation region, 
thus powering its soft X-ray emission. Second, the Lorentz 
torque associated to currents crossing field lines in the pri- 
mary atmosphere (and in the conducting secondary) re- 
distributes angular momentum between the primary spin 
and the orbital motion. In presence of significant GW- 
emission, as expected for two white dwarfs orbiting each 
other at ultrashort periods, there is an additional effect: 
orbital angular momentum is continuously lost from the 
system, causing the orbital period to steadily decrease 
(and the orbit to shrink). Hence, as long as the primary 
is not efficiently kept synchronous by tidal forces, its spin 
will lag behind the orbital motion on the orbital evolu- 
tionary timescale; electric coupling will thus be active and 
exchanging angular momentum between the orbit and the 
primary spin. 

We stress that the absence of any mechanism other than 
the Unipolar Inductor able to affect the primary's spin 
over the orbit evolutionary timescale is an essential as- 
sumption of our analysis: in Appendix A we show that 
tidal synchronization of the primary is indeed not ex- 
pected to be efficient on this timescale, while it may well 
be effective in rapidly synchronizing a low-mass compan- 
ion. The latter effect is discussed in Appendix A as well, 
being however of much lower relevance. 
In this work we consider binary systems consisting of two 
degenerate white dwarfs, with the following mass-radius 
relation (Nauenberg 1972, Marsh & Steeghs 2002): 
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where M is expressed in solar masses. 
In order to make all above features more quantitative we 
need proper expressions for the physical quantities of in- 
terest. Let us start recalling eq. E4 from Appendix E of 
Wu et al. (2002), that describes the evolution of the or- 
bital period of a binary system with the coupled effects of 
UIM and GW: 
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Here Eg = -{32/5)G/c^[q/{l + q)YMla'^ujl is the energy 
loss rate through GW emission (Landau & Lifshitz 1951), 
the second term in parenthesis represents the contribution 
of electric coupling between the primary spin and the or- 
bit and W is the electric current dissipation rate (source 
luminosity). The function g{uJo) < is: 
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The first term (including the coefficient -1/3) corresponds 
to the ratio between Eg and 6jo/^o for two point masses 



with no spin-orbit coupling of any kind. It represents two 
thirds of the total orbital (gravitational plus kinetic) en- 
ergy of the binary system {Eg). The second term in square 
brackets - call it h{a) - is of order unity for most plausi- 
ble system parameters: in particular, 0.6 < h{a) < 1 for 
orbital periods longer than 200 s, primary mass Mi > 0.4 
Mq and secondary mass M2 > 0.08 Mq but, unless ex- 
treme values of all parameters at the same time are as- 
sumed, it is > 0.85 in most cases. This term accounts for 
the secondary spin being tidally locked to the orbit; as the 
system shrinks due to GW emission, an additional tiny 
amount of orbital angular momentum is lost to spin up 
the secondary and the resulting ujo/iOo will be just slightly 
higher, given Mi,q and ujo- Eq- ® can thus be writ- 
ten concisely as g{wo) — (2/3)£'g/i(a), a physically much 
clearer expression which will be frequently used especially 
in § 13 As both Eg and g{ojo) are negative, GW always 
give a positive contribution to lOq while electric coupling 
{W > by definition) may either favour or oppose the 
orbital spin-up depending on the sign of (1 — a). 
Since the two candidate ultrashort period binaries have a 
measured orbital spin-up, in the following we will be inter- 
ested in systems where the orbit is shrinking: if a < 1 this 
is warranted, but when a > 1 spin-orbit coupling transfers 
angular momentum to the orbit and, if sufficiently strong, 
it may even overcome the effect of GW. 
From a general point of view, however, there is no a priori 
reason not to consider systems subject to such a phase of 
orbital spin-down, due to a strong spin-orbit coupling. In 
H3.ll and Appendix B we briefly comment on this situa- 
tion. For the moment, however, we focus only on systems 
where Wo > 0, a condition that, if a > 1, must be im- 
plicitly expressed as \Eg\ > 1 1^/(1 — a)\ and will be made 
more explicit in the next sections. 
The quantity W can be expressed as: 
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where fii is the primary magnetic moment and the sys- 
tem's effective resistance 3? is (see Wu et al. 2002): 
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where N includes G, Mi and q after writing a according to 
Kepler's third law. In the above formula a is the height av- 
eraged WD atmospheric conductivity, H the atmospheric 
depth at which currents cross magnetic field lines and re- 
turn back to the secondary and Ad the thickness of the 
arc-like cross section of the current layer at the primary 
atmosphere. Finally, the geometric factor ]{e) ^ 1 when 
the orbital period is less than 1 hr or so (see Wu et al. 
2002). Combining the two above relations we obtain: 
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where the last equality defines k. 

Finally, combining eq. E4 and E5 from Appendix E of Wu 
et al. (2002), the following expression for the evolution of 
a is obtained: 
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where Ii is the moment of inertia of the primary star. 



2.1. General Considerations and Approximations 

Eq. (UI) in its general form is a non linear, first order dif- 
ferential equation for a, with time-dependent coefficients, 
the non-linear term being due to the coupling between UI 
and GW in the orbital evolution, i.e. by the torque that 
exchanges angular momentum between the primary spin 
and the orbit. 

Independent of the exact solution, a general remark can be 
made that provides a better understanding of the problem 
and allows finding approximations that simplify its math- 
ematical treatment without affecting its very nature. 
The key point is that the orbital period changes contin- 
uously over time: hence, synchronization would be main- 
tained only if the primary spin changed continuously as 
well. If UI was the only mechanism affecting the primary 
spin - as we are assuming - the electrical circuit would 
have to remain always active in order to allow uji to track 

UJo. 

Note that this conclusion holds in general, whatever the 
functional form of lOo is. On the other hand, the details 
of how the system evolves and the state to which it is led 
depend on the functional form of Coo, i.e. on the model 
adopted. 

In the following we obtain an analytic solution to the evo- 
lutionary equation in a specific approximation of the UIM 
and show that it implies the existence of a slightly asyn- 
chronous, asymptotic state^. The physical conditions un- 
der which our particular solution applies are described in 
detail, but we stress that the conclusion concerning the ex- 
istence of a slightly asynchronous, asymptotic state should 
have a more general validity. 

The most natural choice to begin with is to focus on the 
simplified problem in which the effects of UI and GW on 
the evolution of the asynchronism can be decoupled, thus 
neglecting the non-linear term in the evolutionary equa- 
tion. The primary spin evolution is clearly independent 
of GW and is driven by UI alone; so let us specialize to 
the case where GW alone drive the orbital evolution, UI 
providing only a negligible contribution to it. 



The meaning of this will be extensively discussed in the 
rest of the paper. 



2.1.1. Validity of the Approximation 

We check the relative magnitudes of the terms^ in eq. ^ 
to determine the conditions that allow this approximation 
to be introduced consistently. Using eq. (O to re-express 
UJo I UJo in |(7J) it obtains: 
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Let us first compare (B), the contribution of UI to the 
orbital evolution to (C), the contribution of UI to the 
primary spin evolution. The condition _B <C C (say, 
B < IQ-^C) will be met if: 



a < 



10 hiJ^o 



OC UJo 



(9) 



where the right-hand term is a decreasing function of uJo 
and, for component masses in the range (Mi ~ 0.5 1 
M0), gives a < 2 even at an orbital period Po ^ 300 s. 
Therefore, if Porb > 300 s the condition B <C C is met for 
most plausible values of the asynchronism. 
We still need to determine the conditions under which 
B A is true as well and, then, that GW mostly deter- 
mine the orbital evolution. By requiring B < 10^^ A we 
obtain: 
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where c.g.s. units are omitted in the normalizations and a 
height-averaged conductivity = 3 x lO"'^^ e.s.u. has been 
assumed (Wu et al. 2002). Here fii has been normalized 
to a low value because this gives more easily shrinking 
orbits. Indeed, it can be directly checked from the above 
equation that, if /zi ~ 10^^ G cm^ as suggested originally 
by Wu et al. (2002), Pqw > s even for (1 - a) - 0.1. 
This can be stated as follows: relatively highly magnetized 
systems with a > 1 do not shrink unless their degree of 
asynchonism is very low, because spin-orbit coupling dom- 
inates over GW emission. They can shrink only if either 
they are formed with a very close to one (which would 
require some spin-orbit coupling even during the common 
envelope phase that likely leads to their formation) or be- 
gin their evolution with a < 1, quite a remarkable slow 
spin if one considers an initial orbital period ~ lO'* s. 
What is more likely is that at an orbital period of, say, 5 



more precisely, we compare their absolute values 
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hr the primary spin is faster than that (say, a ~ a few). 
In this case, a 10'^^ G cm"^ magnetic moment would cause 

17 /3 

the orbit to widen until (1 — a)ujo is small enough that 
GW dominate over spin-orbit coupling in Wo- 
Let us restrict attention to weakly magnetized systems 
{fii < 10'^^ G cm^): even in this case, condition (|10|l is 
more constraining than ((SJ: it can be met by low mass 
binaries, for fii ^ a few 10'^" G cm^ and P < 2000 s, only 
if (1 — a) is quite small (< 10^^). For higher mass systems 
(Ml > 0.8 Mq), on the other hand, this condition is met 
at short periods even with (1 — a) = (0.1 1). 
The choice of low magnetic moments is completely arbi- 
trary here and justified only because it gives ultrashort 
period systems where lOo is determined by GW alone. In 
paper II we apply the model to RX J0806+15 and RX 
J1914-f 24, showing that such low magnetic moments are 
not simply convenient but seem to be required for the 
model to work, which it does quite well indeed. 
In summary, B can be neglected even for highly asyn- 
chronous systems at orbital periods longer than a few 
thousands seconds and /ii < lO'^^ G cm'^, while at shorter 
periods (100 2000) s this approximation holds only if 
they are almost synchronous (1 — a < 10^^) or have suffi- 
ciently high-mass primaries. In Fig.^we show an example 
of the numerical integration of the evolutionary equation 
Q for two different values of the primary magnetic mo- 
ment. Representative, although arbitrary, system parame- 
ters and initial conditions have been chosen (see captions). 
The early phase of orbital spin-down (the period actually 
changes just slightly) for the highly magnetized system is 
seen clearly in the upper curve of Fig. Q 

2.2. A Steady-State Solution 

At this point we introduce a major simplification of the 
problem by assuming that the timescale over which a 
changes is considerably shorter than the evolutionary 
timescale of ujo- 

On one hand, this restricts the physical regime of inter- 
est to conditions that will be carefully explored in the 
next section. On the other hand it affords a very simple 
and straightforward solution, providing physical insight 
on the problem at hand. An exact, non- linear and time- 
dependent study will be carried out in a future investiga- 
tion (Dall'Osso et al. in preparation). 
With the above approximations, eq. ||7)| becomes a first 
order differential equation for a with constant coefficients 
whose solution is given by: 

a(i) = (ao -aoo)e~*/^° -l-aoo (11) 

where ao is the (arbitrary) initial value and Tq, and 
are defined as follows: 

flllo k ii/3\ 

V^o h ) 
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l-y t — ^OO 
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Time (yrs) 

Fig. 1. Illustrative example of the early orbital evolution 
of a DDB subject to the UI mechanism. The plots are 
obtained through a numerical integration of the coupled 
evolutionary equations for a and ujo- Initial conditions are 
chosen quite arbitrarily, fo illustrative purposes only: the 
initial orbital period is Pi„ = 1.5 x 10^ s (~ 4 hrs) and 
the primary initial spin ~ 25 min {am = 10). Other sys- 
tem parameters are indicated in the figure. Two primary 
magnetic moments were tried and their values, in G cm^, 
label the corresponding curves. 



This solution requires cbo to be driven by GW alone and 
both {ujo/wo) and ojo to be strictly constant: the latter re- 
quires Tq <C To, as both the above quantities evolve on a 
timescale ~ Tq. 

The timescale Tq, and the parameter (1 — aoo) themselves 
are determined by the orbital period and its derivative, so 
they will be subject to secular evolution as well: in par- 
ticular, they both decrease as the system shrinks. Hence, 
once the system has reached the asymptotic asynchronous 
state - starting from an arbitrarily asynchronous configu- 
ration - a will be locked to its "local" steady-state value 
during the subsequent evolution. 

Expressions for Tq, and W can now be rewritten in terms 
of aoo, Ii and the measured quantities tOo and lOq. 

Tq = T^(l - aoo) = To(l - aoo) 
Ulo 

W = aoo/iWoWo-T — (13) 
1 - aoo 

2.3. Validity of the exponential solution 

The picture introduced in the previous section holds un- 
der well defined assumptions that will be addressed here. 
The condition that a changes rapidly with respect to the 

11/3 

timescale over which tOo evolves corresponds to the 
following statement: after a time Tq, luI^^^ must have 
changed by a small amount, say less than 10%, in order 
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for its approximation to a constant coefRcient to be ac- 
ceptable. Eq. H13() then implies: 



We insert in this equation the pure GW expression for 
lOo/lOo and, for the sake of simplicity, write it as Boj^^ , 
where the constant B incorporates all physical constants 
and system constant parameters. It readily obtains: 

BIi , , 

which translates to the following limiting period: 

P < Pf^st = 1.2 X 10^ f ( ) " 

^ \0.9MqJ V2-5xl03oy 

( ^1 y\, ( V {t^) 

where c.g.s. units in the normalizations and j(e) in the de- 
nominator have been omitted. From the above one obtains 
fundamental indications concerning the evolutionary sce- 
nario implied by the UIM: there exists a critical period (in 
the time-dependent case this will essentially correspond to 
a range of periods) shortwards of which the synchroniza- 
tion becomes fast compared to the system orbital evolu- 
tion. Double degenerate binaries are thought to be born 
with orbital periods of a few hours: hence, if asynchronous 
at birth, their asynchronism will not change much faster 
than Ldo until they shrink to a sufficiently small orbital sep- 
aration. All that happens before a system meets the fast 
synchronization requirement depends on its essentially un- 
known initial conditions (initial value of a) and on the 
particular evolutionary path it follows. However, once the 
orbital period is short enough, the UI mechanism becomes 
so efficient as to cause fast synchronization of the primary 
spin, on a timescale Ta <C Tq' in this regime, the value of 
a is brought to the corresponding a^o while oJq remains 
essentially unchanged. 

The state of a system after it has gone through fast syn- 
chronization is completely independent of initial condi- 
tions and its previous evolution; it becomes a function 
of the orbital period and the fundamental parameters 
Ml , q and jii only. 

In particular, higher mass systems reach the fast synchro- 
nization regime at shorter periods than lower mass ones 
(see Fig. because in the latter the current dissipation 
rate is stronger, for a given orbital period and magnetic 
moment and, at the same time, GW are weaker. Hence, 
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Fig. 2. The orbital period for which Tq ^ Tq {Pfast) 
as a function of the primary magnetic moment (/ii) for 
four different values of Mi. The measured uJo and lOq of 
both candidate ultrashort period binaries - RX J0806-I-15 
and RX J1914+24 - constrain secondary masses in the 
range (0.1 0.35) Mq approximately, increasing for de- 
creasing Ml (see paper II): we chose here q in order to 
give an approximately constant M2 ~ 0.2 Mq, for illus- 
trative purposes. The rightmost curve has a significantly 
lower secondary mass, reflecting a similar result of paper 
II, although it would have changed very little if we had al- 
lowed for a somewhat larger M2. A general trend is clear: 
if /ii > 4 X 10^° G cm'^, the condition for fast synchroniza- 
tion is reached for fairly long orbital periods, whatever 
the system masses. Further, the value of Pfast is sensi- 
tive to the primary mass, low-mass systems reaching fast 
synchronization at longer periods than high mass ones for 
this particular choice (and most plausible choices) of the 
mass ratios. Finally, the measured orbital periods of the 
candidate ultrashort period binaries are indicated by the 
dashed lines: if their measured orbital spin-up is due to 
GW alone, they should both be in the fast synchroniza- 
tion regime (or close to it), unless /xi < 10^° G cm"^. This 
is discussed in paper II and is reported here just for the 
sake of illustration. 

high mass sytems can reach shorter periods still maintain- 
ing a relatively high degree of asynchronism, whose exact 
value depends on the initial one. 

Upon inverting relation H17|) (or directly from the example 
of Fig.lSJ it is also found that, if ^1 > 10'^^ G cm^, the fast 
synchronization regime is reached at periods longer than 
lO'^ 10** s, whatever the component masses. Given the 
strong efficiency of UI in such highly magnetic systems, 
they will be characterized by very low steady state values 
of the asynchronism and very short timescales to reach it. 
Therefore, somewhat contrary to intuition, the luminosity 
of highly magnetized systems at short orbital periods will 
be quite low because of the tiny degree of asynchronism 
they can sustain. 

Finally, recall that we have found in § 12.21 that, at pe- 



that can be translated to the following condition on the 
orbital period: 

_ (36/i)* [{H/Ad)j{e)]^ [GMi{l + q)]i 

UJo > UJfast — TT^ — 5 

(2ai?3i?3/2)ii 



V f — 

Ml/ V^c 



(15) 
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riods less than ^ 2000 s, UI is negligible with respect 
to GW in low-mass sytems only if their asynchronism is 
quite low. High-mass sytems, on the other hand, may ful- 
fil that requirement even with a higher degree of asyn- 
chronism. However, from eq. I|17(l it follows that low-mass 
systems reach steady-state at orbital periods significantly 
longer than 10^ s, so they will certainly have a low value of 
(1 — a) shortwards of that. On the contrary, systems with 
sufficiently high-mass primaries (and low /ii) can reach 
periods ~ 10"^ s or less without experiencing fast synchro- 
nization and are ultimately more likely not to fulfil the 
requirement UI << GW at short orbital periods. In Fig. 
I^lwe show the late evolution of the asynchronism parame- 
ter as obtained through a numerical integration of eq. Q . 
The same system parameters of Fig. ^ have been used. 
The behaviour of a in the plot confirms that the system 
never becomes exactly synchronous as it shrinks; the value 
of ckoo has a marked dependence on /ii, in the same way 
expressed by eq. ((T^ . 



1000 

Orbital Period (s) 



Fig. 3. Late evolution of the asynchronism parameter a, 
obtained from the numerical integration of eq. Q, for a 
Double Degenerate Binary with Mi = 1 Mq, q=0.2, and 
two different values of fii . These are indicated in the figure 
and expressed in G cm'^ . The initial spin period is ~ 4 hrs 
and = 10 but, once a = Uoo, its value is independent 
of the previous history of the system: it is determined 
only by system parameters, given the orbital period. The 
dependence of aoo on Uo can be understood here only at a 
qualitative level; in Sethis we obtain it in full generality. 



3. Energy Budget in the UIIVI: powering the 
circuit through GW emission 

There is a general line of reasoning, based on energetic ar- 
guments, that leads to the definition of a system's steady- 
state without solving the evolutionary equation and with- 
out simplifying assumptions'*. 



In the following we neglect the small term h(a) in g^uJo) 
(S|5]l. This term accounts for the secondary spin being tidally 
locked to the orbit and can be easily added to all formulas. It 
is neglected only for the sake of clarity, as is it does not alter 



In the UIM a simple electric circuit is devised, with the 
secondary star acting as a generator and the primary as 
a resistance. The e.m.f. driving the current is due to the 
asynchonism between the primary spin and the orbital 
motion: hence, we can write the total energy available to 
the circuit simply as: 

UIM = - Lol) = E^'ia^ - 1) (18) 



where iJ^^""^ = (l/2)IiUJg is the primary spin energy for 
perfectly synchronous rotation. 

With this definition Ejjim > if wi > uJq while it is 
negative in the opposite situation. However, what counts 
is the absolute value of Euim, so that a more negative 
value corresponds to a larger energy reservoir. 
From the above equation, the rate of change of the total 
available energy is readily obtained: 



Euim = hujoi^ia - tOo) 



(19) 



and we see that when a = I the available electric energy 
goes to zero but has a nonzero, negative first derivative; 
indeed, in this situation lui — because the circuit is 
switched off, while GW still contribute to ojo- Because of 
this, a becomes smaller than 1 and the negative Euim - 
that with a > 1 implied a decreasing energy reservoir - 
now causes the circuit energy to increase again, in abso- 
lute value: GW start feeding energy to the circuit from 
this point on. 

In the UIM, dissipation of the electric current works 
against the e.m.f. and acts to synchronize the primary 
spin with the orbital motion. This is done through the 
torque iVi ~ hiOi, that transfers angular momentum be- 
tween the primary and the orbit. The torque is readily 
obtained remembering that the primary spin derivative is 
(see eq. Ej): 

W 

(20) 



uji a{l — a)IiLj1 

The rate of work done by the torque on the primary - the 
rate of spin energy change £^^p-„ - is just Niuji so that, 
finally, we have: 

(1) _ 



E 



spin 



1 



-w 



(21) 



Since no net angular momentum loss arises from the cou- 
pling, the same torque with an opposite sign acts on the or- 
bit as well, its rate of work being now -E^^/^^'' = —NilOo = 
-E%/a = -W/{l-a). 

Eq. (pn|l and Q can now be subsituted into eq. ((T^ and, 
remembering eq. (|21() . the following is obtained: 



Euim = E 



spin 



o£l p ^UUIM)E^ 



(22) 



This last expression shows that the rate of change of Ejjim 
receives three different contibutions: 



the physical content of the discussion, leaving all emphasis on 
physically more relevant features and somewhat simplifying the 
notation. 
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— the first term, i^sp-„, represents the rate of change 
of the primary spin because of the work done by the 
Lorentz torque. Obviously, this term always causes the 
absolute value of Ejjjm to decrease, since Ejjjm and 
^spin always have opposite signs. 

— the second term is always negative and represents the 
effect of GW emission. Having a negative sign it causes 
the absolute value of Ejjjm to decrease if wi > Uo, but 
in the opposite case it causes to increase and, 
thus, represents a mechanism for injecting energy in 
the circuit. 

This is the only term that would be present if there 
was no spin-orbit coupling at all. Its effect is better 
understood in this framework, indeed. Suppose uji is 
a constant, greater than ujo- since GW cause uJq to 
increase, the ratio a will decrease, thus consuming the 
available electric energy. From the point when a = 
1 on, the further increase of ujo makes a to become 
increasingly smaller than 1, thus powering the circuit 
battery once again: while the loss of orbital energy 
consumes the electric energy when q > 1, it powers 
the electric circuit when a < 1, thus showing that the 
sign of EjjiM in our definition (eq. I18|l has a direct 
physical meaning. 

— the last term describes spin-orbit coupling itself. It 
represents the rate of work done on the orbit by the 
torque iVi, times the ratio between the primary and 
orbital moments of inertia, where the latter is defined 
as lo = Mia^q/{1 + q) and Eg = -{1/2)IoUjI. 

This term always acts to increase Ejjim in absolute 
value, as it can be checked (remember that ii^g < by 
definition) . 

3.1. Conserving the Total Electric Energy 

The above considerations lead naturally to the following 
question: given that only one term out of three in eq. 

always dissipates electric energy, can a condition be 
reached where dissipation is balanced by the terms feed- 
ing energy to the circuit? If this was the case, one would 
expect the flow of currents never to stop and the X-ray 
emission associated to their dissipation never to fade away. 
To answer this point consider again Ejjim in the form of 
eq. (|19|) : we see that the energy first derivative is zero if 
uJia — LVo- Further, we can write the second derivative: 



the orbital motion (a > 0) this condition exists only for 
a < 1, because the primary spins down (wi < 0) if a > 1. 
When the primary spin is antialigned (a < 0), eq. ipn|l 
shows that ui is always positive and again the condition 
Ejjim — cannot be met. 

Therefore, the energy steady-state in the presence of or- 
bital spin-up exists only with the primary spin aligned 
with and somewhat slower than the orbital motion. This 
is not surprising, since we have seen that GW feed energy 
to the circuit only when a < 1. 

When a > 1 and uJo > 0, the energy first derivative re- 
ceives a negative contribution from both terms in eq. H19|l 
and Ejjim is an ever decreasing quantity. A steady-state 
solution with a > 1 can in principle exist only if Uo < 0, 
i.e. when spin-orbit coupling dominates over GW emis- 
sion. If sufficiently strong, this effect could overcome all 
others and allow Ejjim to stay constant or even increase 
while the orbit widens. 

However, a™ = |(i;o|/|(i;i| > 1 must be compared to 
the definition of the two derivatives (eq. El and eq. 120(1 . 
Consider the largest possible value of Wo, the one obtained 
if GWs are completely negligible; even in this case, in 
order for ujo > to hold we must have a < 3(/i//o). 
The right-hand side of this disequality is very hardly 
larger than 1 (one needs quite small component masses - 
Ml < O.SMq), g < 0.2 - and an orbital period shorter than 
400 s. We conclude that steady-state during spin-down is 
very unlikley to exist in real systems. Furthermore, this 
case is of no relevance for the two ultrashort period bina- 
ries RX J0806-f 15 and RX J1914-I-24 and we wiU neglect it 
from here on. For the sake of completeness, in Appendix B 
we show that orbital spin-down in the UIM is necessarily a 
transient phase, after which every system must eventually 
start shrinking. 

3.2. The 'Energy Steady-State" with orbital spin-up 

We can now write uJi and uJo explicitly in eq. ((19(1 taking 
their expressions given, respectively, in eq. (UJ and eq. 
From this we find a simple expression for the value of a 
at which Ejjim stays constant, that we indicate as a™ to 
distinguish it from the steady-state aoo of the approximate 
solution (ESI- 

By imposing Ejjim = it obtains, after a little algebra^ 



EuiM = /iWo(wiQ! - UJo) + IiLOoiioia -\- LUia - Cbo) (23) 

and see that this is zero as well, if Cjia = lOq- Hence, if this 
equality is verified, it implies that the available energy in 
the generator stays constant following the system evolu- 
tion: the dissipation of electric current is exactly balanced 
by the spin-orbit coupling and GW emission. We call this 
the "energy steady-state" . 

We remember here that in the present work we are es- 
sentially concerned with the case when lOq > 0. In this 
regime, it must be Coia > in order for aiui — lUq to be 
possible. When the primary spin has the same verse as 



in ^enN h (^o/^o) . 
.(l-«oo) = ■ 11/3 

kujo 



(24) 



^ We leave {ljo/ijJo) because this is a measured quantity, inde- 
pendent of assumed system parameters. Writing ujo explicitly. 



B 

Then a. 



is obtained as a function of system parameters only. Put 
{'i/2Eg)hLol and C = {'i/2){EJ Eg)h/ {ktol^"'). 

(1/2) 



AC 



Requiring the 



(1-B)± 7(1 + 5) 
expression under square root to be positive, a limit very simi- 
lar to eq. 1281 is obtained, coincident with it if Ii/Io 0, an 
appropriate limit indeed. 



S.DairOsso et al.: Astrophysical UI powered by GW emission 



9 



a second order equation with the foUowing roots: 



1± VT^^ 



(25) 



From the above we see that, in order for a steady-state 
solution to exist, it must he A < 1/4. If A is greater than 
that, the discriminant of the above equation is negative, 
which ultimately implies Euim < 0. 

We address the meaning of the requirement A < 1/4 start- 
ing from the definition: 



1 kujo 

< 

- 4 /i 



11/3 



(26) 



that, remembering the general expression for u!o/(jJo 
and writing the pure GW contribution to it as BwV^ , 
leads to the following: 



be close to 1 (although always somewhat smaller) and the 
other tends to zero (although always somewhat higher). 
Therefore, it seems natural to neglect the latter if the sys- 
tem evolves from a > 1 towards steady-state. If the system 
started from an almost zero initial spin, on the other hand, 
the solution with the smallest a™ would be met first and 
should be taken into account. 

It can be easily realized that the approximate expression 
for aoo (eq. 112(1 differs from the general one for a™ (eq. 
I25|l just by a small quantity. Indeed, from eq. 1(12(1 we can 
write: 



A 



l + A 



(29) 



from which the fast synchronization condition (eq. 114(1 be- 
comes: 



1 - a < 



2EgB 

3kuj^ 6 hujo 



(!-«)< 



1 E„ 



1 kuJn 



4 h 



A 3 

< or 



(27) 



which assures that A < 1/4. 

This condition can in principle be met for any value of a, 
greater or smaller than 1, depending on system parame- 
ters and initial conditions, because it concerns only the 
existence of a steady-state configuration but not the fact 
that the system has actually reached it. However, the con- 
dition must clearly hold when the system actually reaches 
steady-state, and thus has a = a™ < 1: this implies 
n{ujo) > 0, since 1 — is positive, which translates to: 



Wo > 4 



BIl _ UJfast 



(28) 



This corresponds to a limiting period Psteady — ^Pfast 
(eq. 117(1 longwards of which the condition A < 1/4 for 
the existence of an "energy steady-state" is incompatible 
with a™ < 1. Stated differently, this result implies than 
no steady state exists at orbital periods longer than 
Psteady and EjjiM cau only be negative in these cases. 



3.3. Asynchronous steady-state: the general expression 
and its approximation 

We comment here on the relation between the "energy 
steady-state" , obtained in full generality, and the one de- 
rived in the previous section with neglect of the spin-orbit 
coupling and temporal variations of lOq. 
First of all, the "energy steady-state" has two different 
roots while in S 12. 21 we only found one. This is due to the 
non-linearity of the problem: the complete evolutionary 
equation for a has a second order term, which leads to 
the two roots of eq. 1(25(1 . On the other hand, in 12.21 we 
explicitly reduced the equation to the linear form, thus 
excluding one of the two solutions. Concerning this point 
we can note that, once A < 1/4, one solution tends to 



l + A - 110 
3 



A < 



107 



(30) 



In this case, since from cq. (I12II aoo = {1 + A)^^ and A 
must be this small, we can re- write it, at first order in A, 
as aoo — 1 — A. Expanding eq. ((25() at first order in A as 
well, the same expression is obtained, thus proving that 

and aoo are coincident, at lowest order in A. 
Further, eq. I(30() shows that an upper limit on A is not dif- 
ferent from a constraint on the synchronization timescale, 
so that condition A < 1/4 can be given this more intu- 
itive meaning, resulting simply in a less restrictive state- 
ment than the short synchronization condition (eq. I14() . 
It happens to be less restrictive than that just because 
it describes a general result, based on no approximation. 
A < 3/107 corresponds indeed to a very specific case, 
where temporal variations of ujo are so slow that can be 
neglected in the evolution of a (although they ultimately 
determine the very existence of an asynchronous steady- 
state). The result of this section does not rely on that as- 
sumption: it adds to the previous analysis the conclusion 
that a well defined steady-state exists even when temporal 
variations of ujo are non-negligible. Nevertheless, the de- 
tails of if and how the system evolves towards it cannot be 
addressed in the time-independent, linear approximation, 
in the general case. Hence, steady-state is already defined 
for systems meeting the requirement ((28(1 . but these are 
evolving towards it on a timescale that is not much shorter 
than Tq: the fast synchronization regime ensues only at the 
even shorter orbital period Pfast- 

The discussion of this section shows that the steady-state 
degree of asynchronism towards which the system evolves 
corresponds to a stationary state for the available electric 
energy: the system adjusts its parameters as to dissipate 
only the energy that is fed by GW emission and spin- 
orbit coupling, maintaining a somewhat negative, constant 
electric energy reservoir. The latter can be expressed as 
Ek = ""(1 - A + a™) > ""(2 - A). 
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4. Conclusions 

In the present work we have discussed some important 
implications of the Unipolar Inductor Model applied to 
ultrashort period DDBs, that had been overlooked in pre- 
vious works. In particular, we have focussed attention on 
systems with orbital spin-up. The main result of our study 
can be summarized as follows: in the framework of the 
UIM, and in systems whose orbit is shrinking, the dissi- 
pation of the e.m.f. and associated currents can be balanced 
by spin- orbit coupling and mainly, when the primary spin 
is slower than the orbital motion, by the emission of grav- 
itational waves. 

When the asynchronism parameter a is smaller than unity 
GW feed the the circuit battery - at the rate at which the 
orbit shrinks - by driving the primary spin out of syn- 
chronism and through the orbital spin up itself. This is a 
remarkable circuit, in which gravitational energy is "con- 
verted" to electric energy, powering a continuous flow of 
currents. An equilibrium with the energy dissipation pro- 
cess can thus be reached, such that the electric circuit is 
not expected to switch off. 

In the early evolutionary stages - long orbital periods, say 
a few hours - it is likely that the primary spin be faster 
than the very slow orbital motion: the degree of asynchro- 
nism decreases at a rate comparable to that at which lOq 
itself evolves, because the Lorentz torque is initially much 
weaker than the GW torque. 

When the orbital period becomes sufficiently short (P < 
Psteady), ou thc othcr hand, a steady-state solution ex- 
ists independent of the current value of a and such sys- 
tems start evolving towards it. When they reach the even 
shorter period Pfast they enter the fast synchronization 
regime, where steady-state is achieved in a very short time 
compared to Tq. The existence of this regime is due to the 
fact that the rate of dissipation of currents is a stronger 
function of lOq than GW are: hence, a point is reached 
where most of the residual electric energy is consumed 
over a short time - during which GW affect it only very 
slightly - at the expense of the primary spin. From this 
point on the circuit is forced to work in an almost syn- 
chronous state, in which the associated dissipation rate is 
balanced by the energy being fed by GW and spin-orbit 
coupling. 

Therefore, one may say that the electric energy reservoir 
of systems born with a > 1 (and sufficiently weak pri- 
mary magnetic moment to allow the orbit to shrink) must 
initially decrease. They dissipate this energy at a rate de- 
termined by the orbital parameters, magnetic moment and 
degree of asynchronism. Their lifetime is thus a dependent 
variable, fixed by the ratio between the initial Euim and 
the dissipation rate Euim- When they reach a ~ 1 the 
initial reservoir is completely consumed. 
At this point a change in the nature of the circuit oc- 
curs: GW revive it by substracting further orbital energy 
(GW "inject" negative energy in the circuit) and ujo con- 
tinues increasing, thus becoming larger than oji. As this 
happens, the primary spin starts tracking the orbital spin- 



up because of the spin-orbit coupling, although remaining 
somewhat lower than lOo- Since a < 1, the steady-state 
condition Coi — (OojoL can be met, if the orbital period 
has become sufficiently short (eq. I28II . Once steady-state 
is reached, uj\ continues increasing somewhat more rapidly 
than iOo. a becomes a slowly incresing parameter, ap- 
proaching 1 from below over a timescale much longer than 
To itself. However, it cannot reach exactly unity since the 
total energy in the circuit remains constant and cannot 
vanish. 

The system lifetime is now virtually infinite (apart from 
the possible onset of mass transfer) while the degree of 
asynchronism becomes a dependent variable, being fixed 
by the condition that the dissipation rate always match 
the rate of work done by GW and spin-orbit coupling. 
Several details of the UIM deserve further investigation in 
order to better assess its predictions; a number of obser- 
vational constraints that were not considered here should 
still be taken into account. In particular, the simple UIM 
we have appealed to has been debated since its early pro- 
posal for the Jupiter-Io system. Assuming a perfectly field- 
aligned current system implies complete neglect of plasma 
inertia effects. In fact, it has been argued by several au- 
thors (Drell et al. 1965, Neubauer 1980, Neubauer 1998, 
Russell & Huddleston 2000, Saur et al. 2004, Lopes & 
Williams 2005) that "field-aligned" currents could be as- 
sociated to standing (in lo's frame) Alfven disturbances; 
these distort significantly the dipole (unperturbed) field 
lines. The very presence of magnetic interaction between 
Jupiter and its satellite, and the presence of a large-scale 
current flow between them, is not questioned. Therefore, 
overall system energetics, the intensity of the current flow 
and exchange of angular momentum through the currents 
themselves remain essentially valid. On the other hand, 
the way the coupling works may be more complicated 
than assumed here. Plasma inertia effects could, even sig- 
nificantly in principle, affect the efficiency of the coupling 
and the geometry of the current flow, as well as its stabil- 
ity. 

The main goal of the present work was that of demon- 
strating the viability of the model, by showing that the 
UI phase is not short-lived because the emission of GW 
can keep the circuit working at any time. To this aim, we 
have referred to the basic UIM since this appeared to give 
the best physical insight into this problem. Account for 
MHD effects, which may well be of relevance to a more 
detailed description of system properties, is delayed to fu- 
ture work. 

5. Appendix A 

Here we check that the primary tidal synchronization 
timescale is expected to be longer than the GW evolu- 
tionary timescale of the binary system, in the weak vis- 
cosity approximation first introduced by Darwin (1879). 
In this approximation it is assumed that the star shape 
has a quadrupole distortion (to leading order) induced by 
the tidal influence of the companion. In a frame corotating 
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with the star under study, the asynchronism between spin 
and orbit reflects in an apparent rotation of the compan- 
ion at the beat frequency ujb — ui — lUo- If stellar matter 
had no viscosity at all, its shape would instantaneously 
adjust to the present position of the companion and the 
tidal bulge would be aligned with the line of the centres. 
A small but finite viscosity introduces a (small) lag in the 
reaction of the tidal bulge to the changing position of the 
companion, so the tides will lag behind the line of the cen- 
tres by a small time lapse t which measures the viscosity 
itself. The small viscosity approximation translates in the 
requirement that r be sufficiently small, so that the stellar 
shape at time t will be adjusted as to align the axis of its 
tidal bulge with the position the companion had at time 

t-T. 

The expression for the evolution of a component spin due 
to the tidal interaction® with a companion in this frame- 
work was derived by Hut (1981) in Appendix A and is: 



i.(ii) .3i4(5 



1 — a 



(31) 



where in that equation we have assumed a circular or- 
bit (e = 0) with orbital and spin axes aligned (i = 0). 
Here k is the apsidal motion constant of the star (~ 0.12 
for a white dwarf, Verbunt and Hut 1983), is the star 
gyration radius, defined as / = rgME? and T is a charac- 
teristic timescale of the tide, related to r simply by (Hut 
1981): 



T 



GMt 



(32) 



Note the strong dependence of Tt on the ratio {R/a) that 
measures how small is the star with respect to its Roche 
lobe, an intuitive result indeed. Signiflcantly different val- 
ues of Tf are thus expected for components with signifi- 
cantly different masses. 

The lag time r can be related to the mean viscosity of the 
star through (Alexander 1973): 



75 GMl 
2247r i?4 



(33) 



In order to neglect tidal interactions in the evolution of 
the primary spin we must require that Tt be considerably 
longer than the orbital evolutionary timescale (driven by 
GW). Hence: 



^k_if_ fRy \l-a\ ^ 96 (GMl) 



Tr- \ a 



5/3 8/3 
qOJo 



50 



(H-g)i/3 



(34) 



where the right-hand expression corresponds to one tenth 

of Wo/wo expected from GW emission. 

With the above expression we finally obtain (rescaling 



® we consider only the equilibrium tide, whose effects are 
expected in general to be much stronger than the dynamical 
tide (Zahn 1977) 



the orbital period to the shortest known, that of RX 
J0806-hl5): 



-p< 10^ 



321 / 



4/3 



1 



(35) 



where P is expressed in seconds and i?9 is the star radius 
in units of 10^ cm. We need to introduce numbers in or- 
der to check the meaning of this condition: we consider 
the case of a relatively massive (say Mi > 0.8 Mq) pri- 
mary star and a significantly less massive secondary (say, 
M2 < 0.25 Mq), in order to describe the more likley situ- 
ation for the two candidate ultrashort period binaries (see 
paper II). From eq. (|35(l it is seen that the most favourable 
case for tidal synchronization to be efficient is that of short 
period systems, so we consider only them here. 
With a primary mass Mi = 0.8 Mq it obtains, at P '--^ 
1000 s and a — 2 (an extremely asynchronous system 
at such a short period), /I > 10"'^^ g cm""'^ s""'^ in order for 
Tt < IOto, or 71 = 10^^ for the two times to be comparable. 
In a more likely situation, such as the 321.5 orbital period 
source RX J0806-K15, with Mi = 0.8 Mq and a smaller 
asynchronism, (1 — a) ^ 10~^, an even larger value for 
the mean viscosity obtains, Ji > 10^* g cm~^ s~^, in or- 
der for tidal effects to act on a timescale shorter than 10 
To- Note that even stronger constraints obtain when con- 
sidering systems with a larger orbital separation than RX 
J0806-hl5. 

Hence, tides are likely to always have long a synchro- 
nization timescale with respect to the orbital evolutionary 
timescale, unless JI > (10^^ 10""^^) g cm"-"^ s~^, a hardly 
plausible value. 

Indeed, according to Kopal (1968), plasma (or radiative) 
viscosity in non-degenerate stars is at most ~ (10^ lO'') 
g cm~-^ s~^. In a degenerate star, however, it could well 
be of the same order of magnitude or somewhat stronger 
than turbulent viscosity in normal stars with convective 
envelopes, > 10^° -;- 10^^ g cm~^ s~^ (see also Alexander 
1973). Essentially the same conclusions were reached by 
Iben, Tutukov and Fedorova (1998). They found that a 
mean viscosity much larger than 10^^ g cm^-'^ s^^ would be 
required in order for tidal synchronization to be efficient 
over the GW timescale: they referred to the suggestion of 
Smarr & Blandford (1976) in order for a relevant viscosity 
to be obtained, namely that a magnetic field > 10^ G with 
a proper orientation may strongly enhance a WD viscos- 
ity. In particular, they speculate an extreme value ^ 10"'^^ 
g cm~^ s~^ can be reached, assuming an electrical con- 
ductivity ~ 10^^ e.s.u. and a sufficiently deep layer where 
dissipation of tidal energy takes place. We stress that this 
highly speculative suggestion is the only one made in the 
literature for such a high value white dwarf internal vis- 
cosity. 

We focus now on the secondary component: consider the 
representative case of a secondary with M2 ^ 0.2Mq: 
hence, a plausible value JI ^ (10^" 10^^) g cm^^s ~^ 
suffices to make tidal synchronization at least as fast as 
orbital evolution even for long orbital periods (~ 10"* s). 
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Overall, then, it seems at least plausible that secondary 
stars are efficiently synchronized by tidal forces in these 
systems. However, given that some room for doubt may 
be left on this subject, we want to stress here that the 
condition of tidal synchronization of the secondary enters 
our model mainly through the function h(a), of the order 
unity, defined in §|5| 

Hence, as long as the spin of the secondary does not sign- 
ficantly alter the e.m.f. calculated in § [3 the major con- 
sequence of a free secondary's spin in our model would 
be setting h{a) = 1, not a significant change indeed. The 
conclusion of this section can only be that tidal effects are 
at least extremely unlikely to be important in the evolu- 
tion of the primary spin in the systems considered here. 
Althoguh large uncertainties still exist in the determina- 
tion of white dwarf viscosities, all estimates point to it 
being not sufhciently strong to affect the primary spin 
over a timescale at least comparable to Tq. 
We note however that, even if tides acted on the same 
timescale as GW emission (they would be faster than GW 
only i{JI> 10^*, even in the extremely favourable case of 
highly asynchronous systems at a period < 1000 s) this 
would in any case make UIM the strongest spin-orbit cou- 
pling mechanism once the fast synchronization regime is 
reached. Hence, even in the regime where tides may pos- 
sibly become as fast as the orbital evolution (P < 10^ s), 
UIM would be much more efficient in affecting the pri- 
mary's spin. 

Finally, low-mass secondaries are much more easily syn- 
chronized by tidal interactions and it seems that they 
are fully consistent with having a tidal synchronization 
timescale shorter than Tq even at relatively long orbital 
periods. 



6. Appendix B 

We have shown in the text that the state with > 1 
cannot exist for orbital periods longer than a few hundred 
seconds or for Mi + M2 > 0.6Mq or so. Further, we have 
briefly discussed systems undergoing orbital spin-down, 
assuming that they must eventually stop and reverse their 
orbital evolution. However, this latter conclusion has not 
been demonstrated and we address it here. Namely, we 
show that a maximum orbital separation (period) exists in 
this case, beyond which GWs take over spin-orbit coupling 
and the two stars can only spiral-in. We begin writing the 
condition that the orbit keeps widening, that is Wo/i^o < 
or g~^{ujo)[Eg — W/ (1 — a)] < at any orbital period. Put 



_ 2/3 



9 9 



from which: 



En — B„Lul^^^ and thus obtain: 



(36) 



a > 1 



By adding or subtracting 1 from eq. (|37|l we eventually 
obtain: 



1) > 



2B„ 



kujo 



7/3 



1 + 



1 B„ 



2 kwl 



7/3 



(38) 



that would require an ever increasing Euim, in order for 
spin-down to continue. 

Concerning this point note that, during orbital spin-down^ 
this condition can be met if spin-orbit coupling transfers 
more energy to the orbit than that lost by dissipation of 
currents and by the decrease of ujo- Indeed, when a > 1 
{EjjiM > 0) spin-orbit coupling is the only term with a 
positive contribution to Euim fea. l22|l . Hence, in order for 
Euim to increase, remembering eq. (I22f) and the relation 



-E 



(UIM) 
orb 



(cfr. 0, we obtain: 



a < 3- 



E^ 



\E„ 



\E„ 



rpiUIM) 
^orb 



(39) 



where all terms have been written as to be positive. In 
our hypothesis spin-orbit coupling (-Eq^/*'^^) is stronger 
than GWs (Eg): in particular, in the extreme case that 
GWs were completely negligible, one would obtain from 
eq. (|39|l the largest upper limit on a. However, even the 
largest upper limit is too constraing: indeed, neglecting 
the second term on the righ-hand side of eq. H39|l : 



a < 3— 



(40) 



the orbital momentum of inertia lo being defined in 
As already stated in the text, the right-hand side of this 
equation is very hardly greater than 1, especially for pe- 
riods longer than a few hundred seconds, while a > 1 is 
a necessary condition here. So, during orbital spin-down, 
Euim is not expected to increase as the system evolves. 
Thus write: 



1 cx 



(41) 



where 5 is positive and can be smaller or larger than 2. 
The evolution of the asynchronism is thus not fast enough 
for to remain negative. In fact, conditions (|41|l and H38|l 
are not compatible in general: their right-hand sides de- 
fine two functions, yioJo) and w{ijJo), with w more strongly 
dependent on ujo than y. It drops to zero more quickly for 
large Wo and diverges more quickly as lOq becomes small. 
Then, suppose to have a system at a sufficiently short 
orbital period and high a > 1 that spin-orbit coupling 
drives the two component stars apart. As the orbit widens, 



(a^ — 1) cx tj( 



-(2-5) 



, while the condition for orbital spin- 



down to continue requires a stronger, negative, depen- 
dance. A point is reached, at which the degree of asyn- 
chronism is not sufhcient anymore to sustain spin-orbit 
coupling against GWs and the orbit must start shrink- 
ing. Therefore, binary systems can only spin down via this 
mechanism for a finite time, after which they must shrink 



(37) and eventually reach a state with < 1. 
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